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Abstract 

The dispersion relations of collective oscillations of the magnetic moment of magnetic dots ar- 
ranged in square-planar arrays and having magnetic moments perpendicular to the array plane 
are calculated. The presence of the external magnetic field perpendicular to the plane of array, as 
well as the uniaxial anisotropy for single dot are taken into account. The ferromagnetic state with 
all the magnetic moments parallel, and chessboard antiferromagnetic state are considered. The 
dispersion relation yields information about the stability of different states of the array. There is 
a critical magnetic field below which the ferromagnetic state is unstable. The antiferromagnetic 
state is stable for small enough magnetic fields. The dispersion relation is non-analytic as the 
value of the wave vector approaches zero. Non-trivial Van Hove anomalies are also found for both 
ferromagnetic and antiferromagnetic states. 

PACS numbers: 75.75.+a, 75.30.Ds, 75.10.Hk 
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I. INTRODUCTION 



During the last decade the most impressive achievements in magnetism were related to 
fabrication, investigation, and application of artificial magnetic materials (see Ref. [l for a 
recent review). The technologies of sputtering and lithography have progressed to the state 
where the manufacture of nanosize, periodic magnetic superlattices of different types is 
feasible. Among them two-dimensional lattices of sub-micron magnetic particles (so-called 
magnetic dots) attract much attention. These magnetic dots, in the form of circular or 
elliptic cylinders, or rectangular prisms, are made of soft magnetic materials such as Co 
and permalloy,-"- or highly anisotropic materials like Dy- and FePt,- and the dot array 
lattice is usually designed to be quadratic or rectangular. In the array dots are separated 
from each other so that direct exchange interaction between dots is completely absent. 
Thus the dipolar interaction is the sole source of coupling between dots and determines the 
pattern of dot magnetic moment orientations that constitutes the physical properties of a dot 
array. Owing to the absence of exchange, magnetic dot arrays constitute promising material 
for high-density magnetic storage media. For this purposes, the dense arrays (with the 
period of order of 100-200 nm) of small enough magnetic dots with the magnetic moments 
perpendicular to the array plane are optimal, see Ref. S For small enough dots (the 
diameter < 100 nm) the magnetization inside of a dot is almost uniform, producing the 
total magnetic moment m ^> //#, where /is is Bohr magneton. Therefore, this can be 
a new kind of magnetic material with purely two-dimensional lattice structure and pure 
dipolar coupling between large enough magnetic moments. 

To describe the physical properties of magnets in general it is essential to take into account 
interactions having various origins and different energy scales. For adjacent spins the spin 
exchange interaction is almost invariably the strongest. If this interaction is ferromagnetic 
it causes uniform spin ordering over distances that usually substantially exceeds the atomic 
spacing. The dipole-dipole interaction, though as a rule is considerably weaker than the 
exchange for adjacent spins, but it extends to much longer range. The competition of these 
two interactions produces magnetic domain structure with long-range non-uniformity of 
magnetization for usual magnetic samples. 

Models of magnetic moments with the dipolar interaction have been theoretically studied 
for more than 60 years, and many physical properties, lacking in the spin-exchanged systems, 
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are known for those models.—* 1 ^ Note first the presence of a non-unique ground state with 
non-trivial continuous degeneracy that is typical both for three-dimensional lattices^ and 
for two-dimensional lattices of planar— i^ 3 - and three-dimensional dipoles.— Two-dimensional 
systems of Ising-type dipoles with perpendicular anisotropy demonstrate a complicated be- 
havior when subject to an external magnetic field perpendicular to the array's plane. — At 
zero field, the ground state of a square 2D array is chessboard antiferromagnetic (AFM), and 
this state is stable at low fields H < H\ = 2.3M, where M = mo/ a 3 is characteristic field 
of magnetic dipole interaction, and a is the lattice spacing. The ferromagnetic (FM) (sat- 
urated) state is stable at high fields H > H = 9.03M, with a cascade of phase transitions 
between complicated magnetic structures at Hi < H < H . 

The models^— were discussed originally in regard to the description of real crystalline 
spin systems in which the dipolar interaction is predominant, such as some rare-earth mag- 
nets. However, fabrication and experimental study of magnetic dot arrays provide new 
physical systems for the testing of basic magnetism models. For many materials, such as 
compounds with rare-earth ions, granular magnets, and diluted solid solution of paramag- 
netic ions in nonmagnetic crystals, the magnetic properties seem to be similar to those of 
dot arrays. Nevertheless, magnetic dot lattices exhibit physical properties, which are absent 
in all above-mentioned systems. First, dot arrays in contrast to layered crystals are literally 
two-dimensional. The manifestations of long-ranged magnetic dipole interaction are princi- 
pally different for of layered crystals and truly two-dimensional systems.— Second, the scale 
of dipolar interaction of two spins (2/i^5') 2 /a 3 , where /ie the Bohr magneton, a is the inter- 
spin distance, even with large spins like S = 7/2 does not exceed several Kelvins, whereas 
for dots with volume as small as 10 4 nm 3 the characteristic magnetic moment m > 10 3 fisS, 
and the characteristic energy can be comparable to or even higher than the thermal energy 
at room temperature. Besides, for a compound with high density of rare earth ions though 
the exchange interaction is small it is not completely negligible. 

As both FM and AFM exchange interactions of adjacent spins lead to magnetic states 
essentially different from those caused by dipolar interaction, these systems can not be 
considered as purely dipolar. The large magnetic moment of each particle is typical also for 
granular materials, but the magnetic dot lattices are different from the latter by the high 
spatial regularity. It is worth mentioning a new class of materials - molecular crystals with 
high-spin molecules™£ However, they are three-dimensional, and the spin of these systems is 
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much less than the effective spin for magnetic dots not exceeding 10 or 15. It is also important 
that the size of magneto-active part of a single high-spin molecule is small compared to the 
total size of molecule; therefore the dipole interaction is weak. 

Hence, magnetic dot arrays are specific new magnetic materials with purely two- 
dimensional and quite regular lattice structure and long distance dipolar coupling between 
magnetic moments, which are rather large and manifest at high temperature. From the point 
of view of dynamical properties this implies that for magnetic dot arrays the well-defined 
modes of collective oscillations characterized by definite quasimomentum should exist, while 
this is admittedly not the case for granular magnets or dilute solid solutions of paramagnetic 
ions. The direct measurement of the dependence u(k) can be done by the Brillouin light 
scattering method. In the pioneering experiments^— no indication for a band structure 
due to the periodic arrangement of the dots was found, see also review articles.— 1 ^ However, 
more recently the dispersion effects for dense arrays were clearly observed by the Brillouin 
light scattering method,— 1 ^ and time resolved scanning Kerr microscopy— stimulating the 
development of the theory for these systems. 

In first theoretical articles, the finite systems were investigated, sometimes with large 
enough number of dots N, such as iV ~ 10 3 -10 4 , see Ref. |26| or even smaller systems.— 
The analytical calculations were performed using the Bloch theorem which is applicable to 
infinite arrays. The cases of spherical particles,— ~— cylindrical particles in uniform state^i^ 
and in the vortex state,— were also considered, see for review.—) 

In the present work we considered collective modes for a square lattice of rather small 
dots of nearly ellipsoidal shape in which magnetization can be considered as uniform within 
a single dot in the presence of an external bias magnetic field. Only the simplest mag- 
netic states of the magnetic dot system, namely, FM state and chessboard AFM state are 
considered. The long-range character of interaction of magnetic dots leads to unique prop- 
erties of collective excitations in this system, which are absent in both continuous thin films 
and for dipole coupled spins in the 3D lattice. Especially for the FM state it is clear that 
non-analytic dependence of collective mode frequencies on quasi-momentum k, appears. 
Namely, as k — > the spectrum has a finite gap Uq, but the dispersion law is nonanalytic 
u> — y uq + c\k\. For AFM state, the spectrum consists from two energy bands, which are 
connected at the border of the Brillouin zone of the lattice at zero magnetic field, but these 
two bands are well separated by an energy gap at finite fields. For this state the unusual 
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extremum, which is a saddle point in the center of the Brillouin zone is found. 



II. MODEL DESCRIPTION. 

For uniform magnetization inside of each particle the state of the dot array is described 
by the full magnetic moment of each dot, rhf, \mA = mo placed in the square lattice sites 
/ = a(e x l x + e y l y ), where l x , l y are integers. The system Hamiltonian describing the dipole 
interaction of dots subject to an external magnetic field H and taking into account the 
uniaxial magnetic anisotropy for each dot reads^. 

H = 2^ - -3 
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(1) 



Here V=(l-l')/ 



l-V 



k is the anisotropy constant for a given dot, which is assumed 
to be uniaxial with a easy axis e z , perpendicular to the dot array plane. For dots made of 
soft magnetic material, the anisotropy is associated with the dot shape (k > or k < for 
dots oblate or oblong along e z , respectively). For dots made of highly anisotropic magnetic 
materials,-^ the crystalline anisotropy of the material can provide some contribution to 
k. The state of a single dot can be characterized by perpendicular magnetization (fh z = 
±moe z ). However, in the following only the most symmetric case will be considered where 
the external field H is perpendicular to the array plane. 

Although oscillations of magnetization in the dot can be treated in a purely classical 
manner, it is convenient to employ the operator approach.— For dots with homogeneous 
magnetization it is natural to introduce magnon creation and annihilation operators for the 
total magnetic moment. To find a dispersion law in the linear approximation, it is sufficient 
to use the formulae ,~ 

m 3 j = m - gfi B aia r , 

m i,i = + a r )^/gm fi B /2, (2) 

m 2,r = ^( a r~ a r) vV^o/W 2 ) 

below we will put Lande factor g = 2, which in this approximation are similar to the familiar 
Holstein-Primakoff or Dyson- Maleev representations -2. Here hb is the Bohr magneton, 1, 2, 



3 denote the projections on the set of orts, for example, e x , e y , e z , for m = m e 2 or e x , —e y , 
—e z for in = —m Q e z . 

III. FERROMAGNETIC STATE OF ARRAY. 

For the FM state of the dot array (m = e z for all dots) the same set of orts in Eq. (j2j), 
e! = e x , e% — e y and e*3 = e z should be used. In the quadratic approximation over the 
operators at and a,- the Hamiltonian reads 



M 



H + H a 1 



M 



lar > -^—a^ar, t 

/ ' •) Z— / i£|3 / 



<5^0 



3/i B M 



E 

i 



E 

5^0 



(4 + iS- 







5 




5 





2" J J 



a?a' r ? + h.c. 

l+ao 



(3) 

where 5 = l x e x + l y e y is a dimensionless lattice vector, H a = nm is the anisotropy field, 
M = m /a 3 is the characteristic value defining the dipolar interaction intensity and having 
the same dimension as usual the 3D magnetization. The collective modes are introduced 
via states a k and a\ of definite quasi-momentum k 
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where N is the total number of dots in an array. The collective modes are defined by the 



quadratic Hamiltonian over a k , a k which acquires the standard form 



H = 2/i B M ^ A ka\a k + - (B k a\a}_ k + B* k a k a- k ^j 



(5) 



When the form of coefficients A k and B k are established, the collective excitation en- 
ergy e(k) = hui(k) may be found by means of u — v Bogolyubov transformations (see, for 
example,—) and universally reads 



e(k) = 2fi B MyjAl - \B k \ 2 , u(k) = "(My/ A* - \B k \\ 

where 7 = is the gyromagnetic ratio. The concrete forms of A k and B k are defined by 

the distribution of the magnetic moments within the array. For the case of interest (parallel 
ordering of dot magnetization) one can find 



A k = h + (3-^a(0) + ^ rr[())~-a(k) 



B k = 3cr c (fc), 



(6) 



where h = Hq/M and /3 = H a /M are dimensionless magnetic field and anisotropy constant, 
respectively; and the dipolar sums a{k) and a c {k) appear 

a c (k) = a'(k) + ia"(k) = £ ■ e^. (7) 

Such sums naturally emerge for any problem involving interaction of dipoles ordered in 
a lattice. For the two-dimensional case the sums cr(k) and a c (k) have a series of peculiari- 



ties discussed in Refs. 
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33l . Now we note only that the sum a(k) at k = converges for 



large I faster than in the three-dimensional case and has the finite value cr(0) = 9.03362. 
On the other hand, the representation a(k) through the integral apphcable in the three- 
dimensional case (see Ref. 1351 ) is not feasible since the corresponding two-dimensional in- 
tegral J dxdy / (x 2 + y 2 ) 3 ^ 2 diverges as x, y — > 0. Therefore, the regular presentation of 
the static demagnetization field H m through the magnetization components M of the form 

^3 3 

H m = —An ^2 (ejiVjMj), where iVj are demagnetizing factors, J^iVj = 1, fails in the two- 

i=l i=l 

dimensional case. In particular, Yafet and Georgy demonstrated^ that for the atomic 
monolayer the ^-projection of static demagnetizing field H m , is (3/2)cr(0) in our notation, 
and it does not coincide with the continuum theory result for a thin film, H m = —4ne z M z . 
The numerically obtained value of cr(0) = 9.03362 gives (3/2)cr(0) = 1.0783 • 4n, thus the 
difference is nearly 8%. For k ^ the specific feature of the two-dimensional case manifests 
itself more vividly as the appearance of a non-analytical dependence on k, namely, at k — > 0, 

a(k) = cr(O) - k ■ F(k), a c (k) = ■ G(k), (8) 

where k = \k\, G(k) and F(k) are the real functions, analytic as k — > with the limit values 
F(k) — > 27ra, G(k) — > 27ra/3 as k — > and it is invariant under symmetry transformations 
of the square lattice. Thus in the collective oscillation spectra non-analytic features appear 
which are absent in the case of a three-dimensional ferromagnet. From here it is seen 
that in contrast to <r(k) the complex sum a c (k) is not invariant with respect to symmetry 
transformations of the square lattice, but these transformations change only the phase factor, 



and therefore do not have an effect on the frequency. The frequency contains only 
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k-G(k), 



u(k)= 1 \lHf + H Q + ^ (+) {k) 



tff + tfo + f £(-)(£), (9) 



where we introduced the combinations of dipole sums E(±)(fc) 

E (±) (jfe) =<t(0) -<r(jfe)±3 <T c (jfe) . (10) 

— * 

These sums are plotted in Fig. 1 as a function of k for symmetrical directions of the 
reciprocal lattice. The field, H^ s = H a — (3/2)Mer(0) is the z-projection of the effective 
anisotropy field, comprising an anisotropy field for the single particle and demagnetization 
field of the array as a whole, and a (0) (3m /2a 3 ) = 1.08 • (4ttM). 
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FIG. 1: The dependence of the combinations of the dipole sums E(±)(fe) on the quasimomentum 
for symmetric directions of the square lattice, full symbols depict £(_)(£;) = a(0) — a(k) — 3 a c (k) , 
open circles represent £(+)(£;) = cr(0) — cr(fc) + 3 cr c (fe) . k B is the maximal value of the wave 
vector modulus for a given direction, corresponding to the border of the Brillouin zone, k B = 7r/a 
for fc||(0, 1)|| (1,0) and k B = 7r/aand k B = ir^/2/a for fc||(l, 1). 



A. Dispersion relation and stability conditions. 



It is remarked that, the expression for the frequency can be rewritten in the universal 
form 

2 — * — * — * ~| 

to (k) = uj + WintS (+ )(A;) ujq + w int E(_)(/c)l , (11) 



with the parameters uo = ^{H^ + Ho), Ui Qt = cumLR 2 /8cl 3 , where ujq can be thought of 
as the mode gap frequency, uj- int determines the magnitude of the mode dispersion caused 
by interaction, um = 4:7r~/M s is a characteristic frequency of the material, which is ujm= 30 
GHz for permalloy It is worth noting, the same structure ( TTTj) appears for the collective 
mode of vortex precession for the array of vortex state dots.— To define the parameters in 
( TTTj) we used Ui nt = 7M /2 = 7mo/2a 3 and the value of the magnetic moment mo = nM s LR 2 
for the single dot of cylindrical shape of thickness L and radius R. 

Let us discuss the characteristic features of the dispersion relation represented by Eq. ( TTTj) 
in more detail. For small values of \k\ both cr(0) — cr{k) and cr c {k) are linearly increasing 
functions of \k\ for all directions of k, a(0) — cr(k) — > 27ra|A;|, cr c (k) — > 2na\ k |/3, as explained 
in the Appendix. Because of this the first bracket in Eq. (TTTj) in the long-wave limit acquires 
the form H^ t + 47rakM, k =\k\. Note that the multiplier before k, in contrast with that for 
static demagnetizing field, is exactly the same as for continuum films. On the other hand 
the linear components in k compensate each other in the second bracket as k — > and only 
the quadratic terms in k remain (see Fig. 1). Thus the magnon spectra have a peculiarity 
as k — > 0, u(k) = uq + 2ir r jMa\k\ (see Fig. 2), while the value of coefficient of \k\ is exactly 
the same as for a thin magnetic film of saturation magnetization, M and thickness, a, see 
Ref. 
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Note the interesting feature that the lattice constant of the array, related to an 
in-plane space scale, plays the role of film thickness. 

Since the spectrum's behavior depends on the value of the external magnetic field by a 
simple additive way, the sole parameter, determining the form of the spectrum, is the ratio 
A = tuo/^int- First, this dispersion relation is strongly anisotropic. For all values of \k\ inside 
the Brillouin zone the frequency u(k) is increasing monotonically for k parallel to the (1,1) 
direction, but for k along the (1,0) axis the dependence can be non-monotonic, as seen in 
Fig. 2. The oscillations have u 2 > and the ferromagnetic state is stable for weak enough 
interaction, u Q > 1.172w int . Near the point of instability the dependence cu(k) for k \\ (1,0) 
has a minimum near the boundary of the Brillouin zone, which is present until A = 3.6 which 
is also seen in Fig. 2. In all this region of parameter space 1.172 < uo/u)- mt < 3.6 some more 
extrema and saddle points are present inside the Brillouin zone. For weaker interactions, 
Co> int < Co> /3.6, the dependence of u{k) becomes monotonic inside all of the Brillouin zone. 

Here it is remarked that the stability condition for the dot array against small perturba- 
tions is nothing but the condition of positive definiteness of the function u 2 (k). The values 
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FIG. 2: The dispersion law for a dot array in the FM state along some symmetric directions at 
different A = wo/^int (shown near the curves). Here and below we use the notations T, X and M 
for symmetric points of the Brillouin zone (0,0), (0,1) and (1,1), respectively. 

of the sums E(±)(fc) vanishes as k — > 0, and the stability of the FM state against the linear 
long wave excitations is broken if + H < 0, or H < 3<r(0)/2 — H a . After the replace- 
ment 3<t(0)/2 — > 4tt and M — > M s this criterion coincides with the result obtained from the 

— * 

continuum theory for a thin magnetic film. Nevertheless the dependence u(k) is essentially 
different than for thin films. The minimal value of the function u 2 (k) is attained at the 
boundary of the Brillouin zone for the quasimomentum, k parallel to the (1,0) axis. At this 
point, the maximal value of 

-(l/2)£ ( _)(£) = 0.5[3 a c (k) - a(0) + a(k)} 

is equal to 0.5859 = 4tt • 0.0466. Because this instability condition is first realized for a non- 
small k || (1,0), it is stricter than the one found in the long-wave approximation. Finally, 
the ferromagnetic state of the dot array loses stability at H < H c , where 



n M 

H c = — max 

2 



a c (k) -a(0) + a(k) 



— H ~ 



~ 1.125 • AttM - H a . (12) 



This value of the critical field differs from that obtained by the continuum approximation, 
AttM by more than 12 %. Moreover, the character of instability for the array of dots coupled 
by the dipole interaction is principally different that for the continuous film. For the dot array 
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as H < H c the unstable mode has the maximal increment at k — (ir/a)e x or k = (n /a)e y 
corresponding to the transition from the FM state with parallel magnetic moments for all 
the dots to the chessboard AFM state. For continuous films the instability occurs only for 
long wave excitations and leads to appearance of long-period domain structures.— 

B. Density of states. 

The complicated behavior of the dispersion curves play an important role in the formation 
of Van Hove singularities. The Van-Hove singularities are connected with the extrema 
of the dispersion law of quasiparticles; namely, with minima, maxima and saddle points. 
Any branch of collective excitations has at least one of these extrema within the Brillouin 
zone. For the interaction of a finite number of neighboring spins, the dispersion law is 
described by an analytical function, and the function u(k) can be approximated by parabolic 
functions in the vicinity of an extremum. In this case the van Hove singularities have a 
standard form. In particular, for a two-dimensional case the points of minima and maxima 
of uj(k), where ui = uj min or uj = w max , result in a finite jump of the density of states, 
D (uj) = C • (uj — uj min ) or D (ui) = C ■ (uj max — uj), respectively, where (x) is the 
Heaviside step function, (x) = 1 for x > and (x) = for x < 0. Saddle points with 
uj = uj c leads to logarithmic singularities of the form AD (uj) — C ■ In [uj c / \uj — uj c \], and the 
appearance or disappearance of one more singularity of such form at some value of frequency 
has to be clearly seen. 

In our case the structure can be richer due to the long-range character of the interaction, as 
extrema can correspond to an non-standard behavior of D(uj). First, note, the non-standard 
(linear in k) dispersion near the gap frequency, ujq produce much weaker singularities in the 
density of states; for uj > Uo simple calculations give 

D(u) = C(u - uj ) ■ Q(uj - uj ) (13) 

instead of finite jump. In addition, the number of extrema could be more than three. In our 
case all these possibilities can be realized at various values of the parameter A = uo/uj int . 

The analysis of the spectrum shows that for the FM state the dependence uj = u(k) 
always has a minimum in the point T (k = 0); moreover, in the vicinity of this point 
(uj — uq) oc \k\. As well u(k) always has the standard parabolic maximum at the point M 
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(c) 

FIG. 3: The profile of the density of state function, normalized by the condition J D(uj)duj/u>Q = 1, 
for different values of the parameter A = wo/^int- a ) weak interaction, A = 6 > A c ; b) strong 
interaction, A = 1.5 < A c ; c) special case A = A c , where more strong singularity appears. 

(k = (1, 1)). Therefore, near the upper and lower edges of the frequency band, the character 
of the singularities of the density of states D{oj) is universal, D{u) oc (lu — u ) near the u 
(see Eq. (1131) ). and D(u) has a finite jump near the maximal frequency, uv ax - At all values 
of A, the logarithmic singularity of the form of AD (u>) = C ■ In [u c / \u — u c \] is also present 
(see Fig. 3). 

For large A, corresponding to a weak interaction of particles, the frequency grows with \k\ 
for all directions of k. In this case the situation is standard: saddle points are located at four 
symmetrical points of the type of X (1,0), and only the three aforementioned singularities are 
present in the density of state function, see Fig. 3a. However, for small A, A < A c = 3.6 the 
dependence u)(k) for k\ \ (0, 1) is non-monotonous: a local minimum with a standard parabolic 
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dependence is placed in the points of the type of X (1,0), where k = ±ire x /a or k = ±ne y /a. 
Thus for small A four saddle points are symmetrically located inside the Brillouin zone 
on the directions such as (0,1) with \k\ < ir/a. A local minimum leads to existence of 
an additional jump of the density of states having the form D (oj) = C\ + C2 • © (oj — 0Ji m ) 
inside the frequency band for A < A c , see Fig. 3b, whereas the saddle points produce standard 
logarithmic singularities. 

It is also seen that at the critical value of the parameter A = A c ~ 3.6, where two Van 
Hove singularities emerge at some critical value, oj c , the density of states shows a new type 
of singularity 

which is stronger than the standard logarithmic singularity (compare Fig. 3c and Figs. 3a,b). 
To calculate this dependence, one have to take into account the terms of order (k — k ) 4 . It is 
remarked here that the value of the parameter A = u /oj int depends on the external magnetic 
field and can be changed continuously for the same sample during the experiment, thus, the 
observation of this singularity is possible. On the other hand, the non-standard (linear in 
k) dispersion near the gap frequency uq produce much weaker singularities in the density of 
states. In particular, for oj > oj simple calculations give AD (u) — C • (u> — oj )Q (oj — oj ) 
instead of finite jump. This singularity is clearly seen for all values of A (see Fig. 3). 



IV. CHESSBOARD ANTIFERROMAGNETIC STATE. 

In the previous section we investigated the states of the array in which all dots are in 
the same magnetic state (FM ground state of array). Next consider the AFM structure 
where the transition from FM to AFM states can be controlled by a weak external magnetic 
field. The spectral analysis of the dot array in the AFM state is of much interest because 
this transition from FM to AFM states provides the possibility to tune the properties of 
collective modes in the system. 

The method developed here might be easily generalized for the case of simple AFM 
states, which might be described within the framework of a few sublattices. To do that, 
it is necessary to split the array onto different sublattices with the same dot state in each 
sublattice, and introduce different Bose operators for each sublattice. 
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A. Brillouin zone, quasimomentum and magnon operators. 

Let us consider the simplest chessboard antiferromagnetic state, with the magnetic mo- 
ments rhf = (— \) lx+ly m e z in the ground state. This system can be treated as the embedding 
of two quadratic sublattices with the lattice spacing ay/2, as shown in Fig. 4. Let us denote 
the sites of the first sublattice which have rhg = e z as a8, 8 = y/2(8iEi + 82^2), where 81, 
82 are integers and S\ = (e x — e y )/y/2, 62 = (e x + e y )/V2 are basis vectors for sublattices. 
Using orts e x , e y , the vector 8 can be expressed as 8 = 8 x e x + 8 y e y , where 8 X = 81 + 82, 
8 y — 82 — 8±. The sites for second sublattice which have rh^ = —m e z can be expressed as 
a/2, p = v / 2[(a*i + 1/2)£i + (//2 + 1/2)£ 2 ], or p = fx x e x + (J, y e y , fx x = ^ + // 2 + 1, Hy = A*2 — A*i, 
where fj,i and /i 2 are integers. 




FIG. 4: ) The scheme of the sublattices for the particle array with the chessboard magnetic 
order; the solid and dashed lines connect particles belonging to the first and second sublattices, 
correspondingly; dotted lines connect particles from different sublattices. b) the shape of the first 
Brillouin zone, the solid line restricts the Brillouin zone of the sublattice, the dashed line restricts 
the Brillouin zone for the whole lattice. The correspondence between q and q, introduced below, 
is also presented. 

To describe the magnetic oscillations for two sublattices, it is sufficient to use two Bose 
operators as follows, 

m z,s = m o~ 2fi B a^ , m x j = (a$ + aCj ^m n B , 

m y j = i (at - agj ^m [i B ; m Ztjl = 2n B bib p - m , (15) 
m x fi = (bp + &t) ^m fi B , rriyfi = -i (bl - b^j y/m fi B . 

Then the Hamiltonian takes the form, containing the terms describing the interaction of 
oscillations both within first and second sublattices, as well as between different sublattices. 
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The form of the corresponding terms is common to that for FM state, see Eq. fl3]), and we 
do not present this long expression here. 

The translation symmetry within any sublattice allows us to introduce magnon creation 
and annihilation operators for magnetic oscillations through the Bose operators for each of 
the sublattices 



a 5 



g g 



where, naturally, the quasi-momentum q takes values within the first Brillouin zone of the 
sublattice (see Fig. 4). In the basis of the sublattices the vector q = q\E\ + ^2^*2, where 
1 91,2! < Tt/ay/2, and in the basis of the basic lattice the vector q = q x e x + q y e y , where 
Qx = {Q1+Q2) I '"\/2, q y = q x = (<?2 — <Zi)/v2, \q XtV \ < n/a. After this transform, the Hamiltonian 
acquires the standard form H — Yl Hqi where 



g 



(Aj+ h)aia^+ (Aj-tybtbq— 



2fi B M v "" y ' q q v 9 1 qq 

'a 



(17) 



where the following notation are used 

1 v-^ 1 l^e ia & „ 

Aff=y t^-> — --> -w- + /9. 
q ^ W ^ w 24^ |5| 3 



? 2 |£|fi 



3^ (/i, + ^) 2 e^ l^e"^ 
2^ ^ ,f?= 2^W ( } 



There is a simple connection between the sums over the sublattices and the previously 
introduced sums over the whole lattice a{k) and cr c (k), see Eq. (jSJ). For example, a simple 
geometrical transformation gives 

e iaq5 ^ 



^— ' 1 jci: 



|5 1 3 2 3 / 2 
<5^0 1 1 



0" g 



where we introduced the vector q = \f2{q\e x + q2e y ), which is derived from the vector q 
by rotation by the angle it /A and by stretching by the value a/2 (see Fig. 4). It is easy to 
see that for the q values within the first Brillouin zone of the sublattice, the corresponding 
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values of q are within the first Brillouin zone of the whole lattice. By using the rules of 
transition between the sums over the sublattices and the sum over the whole lattice one can 
obtain 

^ = ( 1 -^) <r(0) -2^(^ +A 

C i = —^"M • D <t= \ ("M + ^2"' («)) ' 



(19) 



Thus, the Hamiltonian coefficients describing small oscillations of the AFM state do not 
include new dipole sums different from those for the FM case and could be expressed through 
the sums by means of cumbersome, but simple geometrical transformations. It is necessary 
to note that the Hamiltonian coefficients and D$ in the terms like a^b-q-, responsible 
for interaction between sublattices, are periodic relative to vectors of the reciprocal lattice 
of the whole system, i.e. they have lower symmetry than the coefficients responsible for 
interaction within the sublattice. However, it can be proofed that in final expressions for 
collective mode frequencies the translation symmetry of the reciprocal lattice of period 
7r ja\p2 is restored. 



B. Dispersion relation. 

For diagonalization of the Hamiltonian f fT7|) one can use the generalized Bogolyubov u — v 
transformation and introduce the creation and annihilation operators of magnons of two 
branches, ct, c^* and dl, d$, and the creation and annihilation operators of different branches 
commute, [cq,dq] = 0, [c$, dU = 0. For normal modes 6$ = —iu}^(q)c q -, dg = —iu^(q)d^, 
where (q) and (q) are the frequencies of magnon modes. The generalized Bogolyubov 
transform can be written as 

6? = e«c? + V?cl ? + + Vgdlg . (20) 

Comparing the equations of motion for the operators eg, dg (for example, eg = 
—iuj(_)(q)Cj) and the operators ag, bg [ihdg = [a$, H]), the system of equations for the 
coefficients is presented as a unitary transformation: 
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( (Aff+h) -C ? -D q 



-{Aj+h) F g - 



CI 



-F* 



Dl 

Q 



-Dt 



(A?-h) -CS 



D 9 



■(Ai-h)J 



(21) 



In accordance with (IzTj) the frequencies of collective modes of oscillations of dot magnetic 
moments, are defined by a rather cumbersome formulae 



l£>l 



/i 2 -F 2 



|C| S 



±2 A //i 2 (A 2 - F 2 ) + \AD + FC\ 2 - [Im(CT>*)] 2 



(22) 



but an analysis of the dispersion relation can be done numerically. The dependence of u (q) 
for the two branches, u/-\ (q) and u/+\ (q), for specific values of a magnetic field are presented 
in Fig. 5. 




FIG. 5: The dispersion law along some symmetric directions for a dot array with a moderate 
anisotropy value (/3 = 5) in the AFM state at different values of the magnetic field h (shown on 
the figure). The notations T, X and M are the same as on Fig. 2, but for the Brillouin zone for 
sublattice. 



The given expression is essentially simplified at symmetrical points of the Brillouin zone 
(in the zone center and on the edges), where the expression 

9 



Im(CD* 



2 7/2- 



Re[a*(q)a c (q) 



becomes zero at these symmetrical point. 
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The spectrum behavior essentially depends on the external magnetic field value. The 
important and quite interesting (see below) case of small fields is determined by the behavior 
at h — 0. Therefore, we begin with the study of asymptotics of the spectrum in the absence 
of an external magnetic field. Because of the long-range nature of the dipole interaction, 
the magnon spectrum has specific peculiarities in the center of the Brillouin zone, therefore 
it is reasonable to start with an analysis of the long wave limit k — > 0. 

— * 

To first order in \k\ the sums take the form (7) and after simple algebra one obtains the 
asymptotic form of the dispersion relation at h — 0, 



This expression describes a linear growth of frequency of the upper branch with increasing 
q = \q\ in the vicinity of the center of Brillouin zone. On the other hand, for the frequency 
of the lower branch w_ (q) the terms linear in q are canceled out, and one can expect that 
for this branch the dependence on the components of the vector q should be parabolic. The 
numerical data presented in Fig. 5 indicate this. However, there is also a sharp anisotropy 
of the dependence o;_ (q); namely, in the vicinity of the point q = the frequency w_ (q) 
increases with a q — \q\, if q is parallel to one of directions (1,0), and w_ (q) decreases if q is 
parallel to (1,1). Note by virtue of forth order symmetry of the point q = such behavior 
is forbidden for the quadratic terms. The anisotropy in the expansion of any analytical 
function uj(q) in the q components can appear only due to invariants of the fourth order, 
like (q^ + q y ) or q 2 x q 2 y . Actually, again a non-analyticity, caused by a slow convergence of the 
dipolar sums, takes place here. Our analysis shows, in the vicinity of the point q = the 
dispersion relation can be approximated by the expression 



with the values of the coefficients, a ~0.654 and /3 ~ 1.8. Thus, in the zero field case a new 
type of singular behavior of u(q), namely, the presence of specific saddle point forth- fold 
symmetry, appears in the lower magnon branch. 

In the case of nonzero external magnetic field in the vicinity of the point q — > one gets 




(23) 



ou(q) ^ u + 



(24) 
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the following expression at \q\ <C h accurate to within first order in \q\ 

u ± (q) = jM^j p + (3^2/4) ( ^2 - l) a (0) 

• ^S- (1/2) (y2-l)a(0) + Tra\q\±iH. (25) 

The Eq. ( 125]) shows that the frequencies of both branches of the spectrum in the vicinity of 
the center of the Brillouin zone increase linearly with \q\ at a nonzero magnetic field. 

C. Stability of AFM state. 

A study of the obtained dispersion relation allows a determination of a stability region 
for the AFM state, i.e. a region of parameters at which u 2 (q) becomes negative. The 
problem of finding the stability region for the AFM state is essentially simplified by the fact 
that the local minima of the lower branch of the magnon spectrum are located in the center 
(9(o,o) — 0) or in the corners (<f(i,i) = ^{si + Sq)/ (a a/2)) of the Brillouin zone of the sublattice, 
see Fig. 5. At these points some coefficients of the Hamiltonian vanish, that makes a simple 
analytical consideration possible. 

According to numerical calculations, in the absence of the external magnetic field and 
at small fields, the minimal value of frequency corresponds to magnons having the maximal 
value of quasimomentum, q(i,i) = tt(ei + £2)/( a v2) = ne x /a. Also at this point w±(l, 1) = 
A(i } x) ± yjDfi^ + h 2 . An analysis of tu 2 (q) in the vicinity of this point gives that at zero 
field the AFM state is stable only for enough large value of anisotropy, j3 > (3 CT , 

Per = -~<7 C (tt/o, 0) - 7 ~ 3v/ ^ x (0) « 2.4532. (26) 

Z o 

If /3 > Per, then with increasing of the field the AFM state loses its stability relative to 
perturbations with q*(i,i) at \h\ > h cr ^\ = ■JA^ 1 ^ — D 2 X ^ , see. Fig. 5. 

The analysis also shows that for \h\ > h cv = \f — Fq the AFM state loses the stabil- 
ity relative to small q perturbations. As it is seen from Fig. 5, both these conditions are 
important at different values of anisotropy. Thus, the AFM state is stable relative to arbi- 
trary small perturbations subject to the condition, h < h CT = min{/i crj o, ^ C r,n}- For small 
anisotropy /3 < fix — 3.358 an instability with maximum values of q is developed, whereas 
at (3 > (3i the AFM state is unstable relative to long-wave perturbations with q \ja. 
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2 3 4 5 (3 



FIG. 6: The dependence of the critical magnetic field on anisotropy. Here both characteristic fields, 
hcr,o an d /i C r,ii are plotted as functions of /3, and the stability region of the AFM state is doubly 
patterned region under both curves. 

D. Density of states. 

Let us consider positions and form of Van Hove singularities for the spectral density of 
magnons of both branches as shown in Fig. 7. In a weak field, frequency bands corresponding 
to the upper and lower branches of collective oscillations overlap. In this case it is convenient 
to introduce the partial density of states as shown in Fig. 7a for zero field. As for the 
aforementioned FM state of the array nonstandard behavior of D[oj) for both branches can 
merge caused by long-range character of interaction. As the previously discussed example, 
when the discontinuity of the derivative of density of states has the form D(u>) = C 1 + 
C2(ou — uo)Q(co — ojq), this corresponds to the linear dependence of ou(q) in the vicinity of 
the center of the Brillouin zone, u(q) — lu oc \q\. For the AFM case at various values of 
magnetic field this possibility can be realized, but the situation can be richer compared to 
the FM case. In particular, this peculiarity is manifested differently in the density of states 
for the upper and lower bands. 

It happens that the standard scheme of extrema (a minimum in the centre of the Brillouin 
zone, maxima on the edges of the Brillouin zone, saddle points in for symmetrical points 
like X(1,0)) take place only for the upper branch and only at large enough field. In this case 
in the upper band there are only three Van Hove singularities (see Fig. 4b, c). For the upper 
branch the point q = at any field always defines the global minimum and peculiarities like 
D(uj) oc (uj — uo)Q(oj — uo) define the density of states near the lower edge of the upper band. 

— * 

For all fields the absolute maximum of the function co + (k) having a parabolic dependence 
is located at the corners of the Brillouin zone at points such as M (1,1). For this reason, 
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0.0 0.2 0.4 2.2 2.4 2.6 2.8 3.0 3.2 

(c) h=5. 



FIG. 7: The density of state for two branches of magnons for the AFM dot array with /3 = 5, 
normalized by the condition f D(oj)doj/co a = 1, co a = jH a , for each band, for different values of the 
magnetic field, a) the field equals zero, pay attention that these two branches are overlapped, b) 
h = 2.5, at moderate field values a linear peculiarity of the density of states is present inside of the 
lower band of collective modes; c) h = 5 presents the case of large fields, where a linear peculiarity 
of the density of states corresponds to lower frequency, and the finite jump is present inside the 
low energy band. The characteristic frequencies are mentioned by vertical arrows. 

one expects the standard peculiarity of D(u) like a finite discontinuity near the maximal 
frequency. For the upper band a nonstandard behavior manifests itself only at small fields 
h < 1.2, when at points like X(1,0) minima appear, while saddle points move along directions 
like (1,0) into the Brillouin zone. The situation here resembles that described for the FM 
case and we do not discuss it here. 

The scheme of Van Hove singularities for the lower branch are much less standard. First 
of all, for the low band at all values of field the global maximum is located at the points like 
X(1,0). A standard peculiarity of the density of states like a finite jump corresponds to this 
frequency. At h ^ (the zero field case is special and will be considered separately) both at 
the edges of the Brillouin zone [points like M(l,l)] and in its center the minima are present, 



21 



and the relative depth of the minima at points T and M is defined by a field value. Thus, all 
symmetrical points of the Brillouin zone, T, X and M are occupied by minima or maxima. 
For this reason the saddle points, to which standard logarithmic peculiarities of density of 
states correspond, are moved inside of the Brillouin zone along a direction like (1,0). 

The relative position of minima is defined by a value of the magnetic field. For the 
lower branch at weak fields (h < 3) the absolute minimum of the function U-(k) is located 
at a point like M. The parabolic dependence oo(q) corresponds to this minimum, and the 
latter defines a standard peculiarity of the density of states like a finite discontinuity (see 
Fig. 7a,b). At such values of the field the peculiarity connected with the point T (q = 0) 
lies within the lower band frequency and manifests itself as a derivative jump D(u) (for 
example, at u — 0.6u) a for h = 2.5) and it is almost invisible in the figures. For higher 
fields (h > 3) the frequency minimum at the point T (q = 0) becomes deeper, i.e. the local 
minimum value of the lower branch of the spectrum (u;_(M)) lies higher than the frequency 
of long- wave oscillations w_(0), (see Fig. 7). In this case a linear behavior of the density of 
states D{u) ~ (u — ojo)Q(uj — Uq) is clearly seen at the lower edge of the frequency band for 
the lower branch. In this case, within the band there is one more singularity of the density 
of states, which is a finite jump of the density of states, determined by a local parabolic 
minimum at the point M (1,1), see Fig. 7c. 

One of the most interesting singularities of the spectrum of the lower branch is observed 
in the special case of zero field. First, in this case the global maxima at the points X (1,0) 
have a non-parabolic dependence u(q), u(q) = u) max — c\Aq\, where Aq is a deviation of q 
from the point X. For this reason, the upper edge of the energy band is characterized be 
non-standard linear Van Hove singularity, D(u) = C(w max — u) ■ 0(w max — uj). On the other 
hand in the expansion of oo(q) in the center of the Brillouin zone the linear term in q is absent, 
and there is a specific saddle point with four-fold symmetry having a non-analytic behavior, 
see Eq. ( |24l) . However the analysis shows, that this saddle point leads to the standard 
logarithmic Van Hove singularity. For small, but finite value of the external field in the 
zone center there is a linear non-analytical behavior of oj(q) ~ ujq + a\q\, which is typical for 
the systems having dipole interaction. While at the point q = there is a minimum and a 
non-analytic saddle point splits into four standard saddle points moved towards the points 
X. 
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V. CONCLUSION. 



First let us discuss one possible generalization of our approach for the problem of interest 
not considered here. The concrete calculations here have been done for the most symmetric 
configurations of the system; namely the ideal square lattice array, and individual dots 
having rotational uniaxial symmetry. The consideration of the same problem with lower 
symmetry will not produce any principal difficulties. For a lattice of lower symmetry, such as 
rectangular, the same peculiarities for dipolar sums and long wave asymptotics of dispersion 
relations are present. The case of dots of lower symmetry, such as rectangular or elliptic, 
with almost homogeneous magnetization within the dots can be easily analyzed, just by 
adding to Eq. ([T]) extra terms describing more complicated (biaxial) anisotropy of a single 
dot. Qualitatively the results will be the same. 

In conclusion, the dynamic properties of the array of magnetic dots of different types 
without direct exchange interaction between dots have been investigated. For such systems, 
only the magnetic dipole interaction can be a source of dot interaction. Exploiting the 
translation symmetry of the array, and by use of the quasiparticle (magnon) formalism, we 
have found full spectra for the change of quasimomentum within all of the Brillouin zone. 
The direct measurement of the dependence u(k) can be done by the Brillouin light scattering 
method. Our calculations demonstrate that the dispersion is a strongly increasing function 
of the dot thickness. It will be interesting to observe non-monotonic dependence u(k) and 
especially, the change of the character of this dependence with the change of the magnetic 
field found here. The most impressive change of the oj(k) dependence should appears near 
the transition from the FM to AFM states of an array. This property can be of interest 
for the design of a new generation of microwave devices in the modern direction of the 
applied physics of magnetic nanoparticles, i.e., the so-called magnonics, which has been 
widely discussed in the literature.— 

An important application of the calculated dispersion laws is the investigation of the 
stability of given magnetic states of the array. It has been found that for arrays with 
homogeneous magnetization within a dot directed perpendicular to the array's plane there 
will always be a nonzero critical field, either external or in combination with an anisotropy 
field. For ferromagnetic state of array, the instability is developed for the quasimomentum at 
the points of the type of (1,0) at the border of the Brillouin zone, and it leads to a transition 
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to the AFM state with chessboard orientation of the dot's magnetic moment. 

For all magnetic states of the array investigated here, the dispersion relation is non- 
analytic as k — > because of the long-range nature of the dipolar interaction of oscillating 
magnetic moments of dots and the presence of singularities in dipolar sums. For such 
modes with finite gap frequency, u the magnon spectra have the peculiarity of the type 
u(k) — uq oc \k\ as k — > 0. For the AFM structure without a magnetic field, a new sort 
of singularity, leading to saddle point with four-fold symmetry of the function u(k) near 
the value k = is predicted. An important and non-trivial property of most of the modes 
considered here is the strong anisotropy of the function u(k). Note also the decreasing or 
non-monotonic dependence of the mode frequency oj(k) on the wave vector k, observed for 
non-small interaction of dots. 

The non-standard dependence u(k) for small k produce the change of the density of 
states D(u). For the spectrum of the FM state, instead of standard two-dimensional Van 
Hove singularities of a form of finite jump near the gap frequency, uq, or singularities like 
In (a> c /|u; c — oj\), weaker singularities of the form D(uj) ~ (u — uq)Q(u — uq) appear here. 
On the other hand, for dispersion relations for the FM state ( iTTj) the change of the character 
of its extremum at the point X of the type of (1,0) is predicted at some ratio of the system 
parameters, uo/ui nt — 3.6 (see Fig. 2). For this critical value of the parameter uo/uJi n t the 
density of states has the singularity D(u) oc l/\cu — cod 1 ^, where co c is the value of frequency 
at this critical point, which is stronger than the standard finite jump. The value of the 
parameter wo/^int depends on the external magnetic field and can be changed continuously 
during the experiment. 
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Appendix. 



The expressions for collective mode frequencies contain series such as the dipole sums 
cr{k) and a c (k). Here and below in this Section we will use the dimensionless vector I and 
the condition I ^ in the sums is implied. The mathematical properties of such series 
are important not only for this problem, but also for any example of lattice systems with 
identical particles coupled by the dipole interaction. 

Let us discuss properties of these series. As we will demonstrate these series have very 
singular behavior as a functions of quasimomentum k, and the double sum converges rather 
slowly. This manifests in sum properties near symmetrical points of the reciprocal lattice k , 
especially near the point of origin k = or near the points of the type of M (1,1), for which 
fco = ±(e x ± e y )Tr/a, and X (1,1), for which ko = ±(e x )n/a or ko = (±e y )-n-/a. Analyzing 
small deviations from these points, k = k + q, where q is small, one normally has to calculate 
derivatives such as [d 2 a(k)/dqidqj] at the point k = k . Term by term differentiation of the 
dipole sums gives series like E^ (Vl^l) ' exp(ifcoZ), which are alternating and converge only 
conditionally. Moreover, for k = 0, i.e., for the physically most interesting case of long 
wave oscillations, the corresponding coefficient of ef 2 is described by the divergent series 
The same property is present for the complex sum a c (k), which is also important 
for the description of dipole-coupled modes. Therefore, the dipole sums at k ~ can be 
non-analytical and it results in a non-standard dispersion relation for oscillations described 
above. 

— * — * — * 

Let us start with study of the sum a{k) for small \k\/ks- To analyze the behavior of cr(k) 
near the point k = one can write k = q, \q\ <^l/a and present this series as 

E 5 = E ■ V f - e ~° p + eii - 1 1 - • < 27 > 
l'1 l'1 

where the multiplier exp(— a\l\ 2 ) is chosen to provide a fast convergence of the corresponding 
series. The first term on the right hand side of Eq. fj27|) converges rapidly for |/| > 1/i/a 
and at q — \q\ — > it contributes as —q 2 D', where D' = (1/|/|) exp (-ah is an absolutely 
converging series. Non-analyticity as q — > is defined by the second term, calculation of 
which can be simplified. Indeed, at |/| <^ 1/y/a it comprises the small multiplier af 2 , and 
the contribution from the region I < 1/y/a is expected to be small. For the outer region 
|/| > 1/y/a one can expect that discreteness effects are small and the sum can be replaced 
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by the integral, Y^,T ~* j rdrdx, where r and x are the polar coordinates. Integration over 
X can be done exactly through the Bessel function Jo(qr) and the non-analytical part of the 

oo 

sum cr(k) at small \k\ is determined by the integral I = 2tt j (dr/r 2 ) J$(qr) (1 — exp(— or 2 )), 

o 

which is convergent both as r — Y and r — > oo. It is convenient to integrate by parts, which 

oo 

in the case we are interested in (a — » ) yields I = 2n j (dr /r)J\(qr) = 2ir\q\. From here 

o 

we arrive at the expression, used in <^ 

a{k) = tr(O) - 27ra\k\ at k -> . 

For the complex sum a c (k) near the origin k — q, \q\ -C 1/a the same approach gives the 
double integral, 

OO 27T 

<7 c (g) = y ^ • y dxe~ 2lx • e igrcos(x - a) , 
o o 

where a is the angle between q and the a? axis. Here the integral over x can be done, which 
gives <j c = 2ire~ 2ia J(dr/r 2 ) ^(gr). Further, for the integral over r, which is convergent 
both as r — )• due to asymptotic ^(z) — ^ (1/2) {z/2) 2 as z — )• and for r — >■ oo, it is not 
necessary to include the regularization multiplier exp(— ar 2 ) like in (1271) . Finally we arrive 
at the following expression 

r 2tt (k x -ik y ) 2 r 
<J c \k) = — a 1 — — as k — > 0. 

Thus the relation between er(&) and a c (k), cr{k) — >■ 3|cr c (fc)| as — )■ 0, which is of great 
importance for the spectral analysis is asymptotically exact at small values of \k\. This 
result is in accordance with the numerical calculation. 

Let us discuss now the calculation of the series cr(k) and a c (k) near the second symmetric 
point k of the type of (1,1). Here for both series the values of = [d 2 a(k) / dkidkj]j: = ^ o are 
described by convergent alternating series like $^P2p(0 ■ (—l) lx+ly (l 2 +ly) < - 2p+1 ^ 2 ) where 
P2p(l) is the polynomial in l x , l y of degree 2p, p = 1 and p = 2 for components Dij in the 
case of cr(fc) and a c (k), respectively. 

By using symmetry relations, one can demonstrate that the value of the function a(k) at 

~~ * — 1 /2 

k of type (1,1) is described only by one sum, D = — (—l) lx+ly (Z 2 , + I 2 ) , the value of 
which can be easily found numerically, D = 1.61554 . As a result the following expression 
is obtained 

a(k) £ a(k Q ) + ^Da 2 \k - k \ 2 . 
28 



For analysis of the sum a c (k) near the point k = k we note that because of the evident 
relations 

the value of a c (ko) = 0. The next terms of the expansion a c (k) over the vector components 
q = k — k can be written easily, 



[ l l + l y) 



They comprise both real terms proportional to q 2 — q 2 and imaginary terms like q x q y . 
After simple algebra, these terms are presented via the sum D introduced above and one 
more convergent sum of similar structure 

D 1 = - - l 2 y ) 2 (-l) lx+ly l (l 2 x + l 2 y )' 5/2 = 3.31891. 

Finally, in the quadratic approximation over q = k — k near the point k of the type (1,1) 
one can present the sum as 

v c (k) = cr c (k ) + i{D 1 - D)a 2 q x q y + ^Dxa 2 (ql - q 2 y ). 

This function, in contrast with a{k) near the same point ko, is anisotropic. However, 
expressions for the frequencies contain a(k) and a c (k) only in the combination (J 2 (k) — 
9|(j c (A;)| 2 , therefore for presentation of the spectrum of collective modes accurate to \q\ 2 it 
is not essential to take a c {k) into account, and to the accuracy of q 2 the spectrum near the 
point k of the type (1,1) is radially symmetric over components of the vector q . 

The point with ko of the type (1,0) including the points like ko = ±7ie x /a or k = ±7ie y /a, 
possess lower (biaxial) symmetry than the points k = or points of the type (1,1). The 

— * — * 

rather complicated study in the vicinity of this point (q = k — k is small) shows the 
analytical dependence for the components of q, in this case for the components of q parallel 

— * 

and perpendicular to ko 

a(k)~a(ko) = ^a 2 (D liq 2 -D ± q 2 ± ), 
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where D\\ and D± are convergent series of the form 

/2(1 )lx 

D \\ = -Y, 5/2 - - 6354 > 



£) ± = \^^ ' ^ 1.256. (29) 

^ C/2 i /2A 5 / 2 v ' 

fyo Vx ^ l y) 
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